Hesitant 2-tuple fuzzy linguistic multi-criteria decision-making method based on correlation measures

Correlation is considered the most important factor in analyzing the data in statistics. It is used to measure the movement of two different variables linearly. The concept of correlation is well-known and used in different fields to measure the association between two variables. The hesitant 2-tuple fuzzy linguistic set (H2FLS) comes out to be valuable in addressing people’s reluctant subjective data. The purpose of this paper is to analyze new correlation measures between H2FLSs and apply them in the decision-making process. First and foremost, the ideas of mean and variance of hesitant 2-tuple fuzzy linguistic elements (H2FLEs) are introduced. Then, a new correlation coefficient between H2FLSs is established. In addition, considering that different H2FLEs may have different criteria weights, the weighted correlation coefficient and ordered weighted correlation coefficient are further investigated. A practical example concerning the detailed procedure of solving problems is exemplified to feature the reasonableness and attainability of the proposed technique in situations where the criteria weights are either known or unknown. When the weight vector is unknown, the best-worst method (BWM) is used to acquire the criteria weights in the context of a hesitant 2-tuple fuzzy linguistic environment. Furthermore, a comparative study is undertaken with current techniques to provide a vision into the design decision-making process. Finally, it is verified that the proposed correlation coefficient between H2FLSs is more satisfactory than the extant ones, and the correlation coefficient with the weights of criteria being either known or unknown is applicable.


Introduction
We mostly interact with multi-criteria group decision-making (MCGDM) problems in our routine life, which are used to choose the best alternatives among different alternatives. Often we deal with dubious situations, and due to the ambiguous nature of MCGDM problems, the fuzzy method may be considered as a valid technique to define the characteristics like linguistic variables [1][2][3], interval values [4,5] and hesitant fuzzy numbers [6]. We have to deal with uncertainty at different levels. The two most important approaches introduced to investigate, interrogate, and tackle this type of uncertainty are fuzzy set theory and probability theory. For the first time, Zadeh [7] characterized the fuzzy linguistic techniques to manage the decisionmaking problems, which utilizes linguistic factors whose contents are not quantities, but phrases or terms, to address subjective assessments. After the development of the linguistic approach, many researchers took an interest in this area and applied it to MCGDM problems [8]. Over time, researchers feel some limitation in the linguistic variable: it can only express the membership degree. There are many situations where consideration of membership and non-membership functions is necessary. For some dynamic issues, it is simple for decision-makers or experts to communicate their assessments in phonetic or linguistic terms, for example, "low" cost, "high" quality, "great" execution. The linguistic variable upgrades the possibility and adaptability of decision models, and the advancement of analyzing linguistic variables have prompted a new research territory named computing with words [7]. However, the customary linguistic methodology uses a oneword linguistic term to evaluate a linguistic variable. This way cannot represent the hesitant impression of decision-makers. For example, while assessing a research proposal, a decisionmaker may state, "it is between good and excellent". Also, in a group decision-making process, the gathering's appraisals could be uncertain if there are various assessments among many individuals. To overthrow the inefficiency of the linguistic process, Rodríguez et al. [9] introduced the concept of hesitant fuzzy linguistic term set (HFLTS), where decision-makers can give their ambiguous evaluations using several linguistic terms. Consequently, the HFLTS has pulled in numerous scholars' attention [10,11].
Herrera and Martins [12,13] proposed a 2-tuple linguistic representation model that can deal with verbal and numerical data viably in decision-making without mutilating and harming data. This helpful model is basically founded on the symbolic translation of linguistic variables and has been generally utilized in different MCDM problems [14]. In light of the difficulty of items and dubiousness of individuals thinking, it is increasingly suitable to give an evaluation by using linguistic data rather than numerical numbers [15,16]. Accordingly, the number of analysts dealt with linguistic data for taking care of MCGDM issues and characterized various techniques to explore the linguistic assessment components. There are mainly three parts of the linguistic information model [16,17], procedure based on membership degree, procedure based on linguistic symbols, and the procedure based in linguistic 2-tuple. However, the initial two techniques have a few restrictions that they played out their interpretation procedure as an estimation procedure to assess the finding in the original space, thus invigorating the absence of exactness and precision [18]. Contrarily, the precision of linguistic terms, calculation, and interpretation of results can be improved with the assistance of the 2-tuple linguistic model. In light of its recognized quality and viability in managing the linguistic judgment of experts, the 2-tuple linguistic models have been generally utilized in decisionmaking problems [19][20][21]. Zhang et al. [22] examined the transitivity of 2-tuple to tackle the issue of group decision-making problems. Apart from this, the idea of a hesitant 2-tuple fuzzy linguistic set is proposed by Wei and Liao [23]. Beg and Rashid [14] extended fuzzy TOPSIS for hesitant 2-tuple linguistic term sets with the assessment of some experts about the properties of alternatives. Wang et al. [24] built up a hesitant 2-tuple linguistic aggregation operator for structuring a new type of group decision-making approach to deal with addressing the data combination, including the interrelationship between aggregated terms and the prioritization relationship among experts under hesitant 2-tuple linguistic circumstance.
At whatever point we talk about the strength of the relationship between two variables, the correlation coefficient is a suitable measurement tool. In all estimation procedures, the correlation coefficient is an effective apparatus and broadly utilized in statistical analysis to quantify the strength or level of relationship between two variables [25][26][27]. However, it is conspicuous that the idea constructed on daily life problems is uncertain. So as to handle this issue, correlation is expanded by different analysts to fuzzy correlation [28][29][30]. Tyagi [31] broadened the idea of the hesitant fuzzy set to dual hesitant fuzzy set and built up a correlation coefficient between dual hesitant fuzzy sets. They applied it to assess the correlation coefficient estimation among various water parameters of four distinct lakes. A correlation coefficient based on the central interval is proposed by Saneifard [27]. Chen et al. [32] outlined a correlation coefficient for hesitant fuzzy sets and applied it to clustering analysis in a hesitant fuzzy environment. Moreover, they likewise proposed some interval-valued hesitant fuzzy sets, introduced their correlation coefficients, and determined their applications in clustering. Arora and Garg [33] took inspiration from the soft set theory and built up the correlation coefficient to measure the relationship between dual hesitant fuzzy soft sets, named the novel correlation coefficient.
MCDM is a developing field of study that helps experts distinguish the best alternative from a set of alternatives assessed through a set of criteria weights [34]. A basic portion of this issue is the arrangement of the significance of criteria weights. Various methods have been developed to recognize these criteria, including [35][36][37], SMARTS and SMARTER, analytic hierarchy process. BWM is also linked to them, designed by Rezaie [38] that depends on an organized pairwise comparison framework for an optimal model to get the ideal weights of criteria. This method has fascinated scholars' attention from different fields since its presentation. Hadi et al. [39] presented z-number BWM to deal with the uncertain situation of MCGDM problems. Mi et al. [40] presented some questions relating to BWM in the form of why what, what for, and what next. Muhammad & Razaei [41] introduced a Bayesian BWM to find the aggregated final criteria weights of models for a group of decision-makers. Furthermore, many other works are available in the literature on BWM, including multiplicative BWM [42], BWM MULTIMOORA [43], BWM VIKOR [44], and others [45].
Adopting inspiration from the existing approaches [24,46], a new type of correlation coefficient between H2FLSs is proposed in this paper. Numerous methodologies exist in the literature that works in information energy and develops a recipe for computing the correlation coefficient. However, there was no technique for estimating the distance-based correlation coefficient from a statistics point of view in the context of H2FLSs. The proposed correlation coefficient can determine how strong the relationship between the two alternatives/variables is. It reveals information about the direction of the relationship and can also display data behaviour that is difficult to study by simple observation. However, drawbacks are also related, such as the inability to determine the curvilinear relationship. Furthermore, it is unclear which alternatives/variables are most affected in the correlation process and which ones are dependent and independent alternatives/variables. It also has no control over the third alternatives/ variables influencing the relationship.

Research challenges and gaps
Correlation coefficients are utilized to show, in a linear design, the closeness level of a couple of variables. The correlation measures can be characterized into two classes, i.e., the statistics based correlation whose findings have a place in the interval [−1, 1] [46], and the information energy-based correlation measures, whose results vary inside the interval [0, 1] [26,30]. Many researchers developed different types of correlation coefficients, as discussed in the introduction, but we still need a statistics-based correlation coefficient that deals H2FLSs. For example, for hesitant fuzzy linguistic term sets, [26,46] suggested correlation coefficients, while [30] proposed a correlation coefficient for dual hesitant fuzzy linguistic term sets. Aside from these methods, there are other alternative correlation coefficients for dealing with various extensions of fuzzy sets, but none was for dealing with H2FLSs. Therefore, there was a challenging work that no technique for estimating the correlation coefficient from a statistics standpoint in the context of H2FLSs is available.

Contributions of this study
In this study, we explore the correlation measure between H2FLSs from a statistical perspective. The theoretical commitments of this paper can be summed up as follows: • We proposed the idea of the mean and variance of the H2FLEs.
• We propose a new correlation coefficient between H2FLSs, whose findings lie in the interval The proposed correlation coefficient between H2FLSs can be utilized to show the positive or negative correlation between two variables in the context of a hesitant 2-tuple fuzzy linguistic environment. It overcomes the limitation of the previous correlation coefficients proposed in [26,30], accordingly is significantly more persuading and has more extensive application possibilities. Practically, the proposed correlation coefficient can be applied in many real-life group decision-making problems, like engineering, medical studies, business studies, etc. The proposed approach is simple to implement, but there are some limitations. For example, the calculations may take longer to complete if the number of alternatives is large during the computation. Furthermore, unlike the conventional correlation coefficient, the proposed correlation coefficient is not independent of origin and scale. It also does not indicate the cause and effect relationship between the alternatives.
Furthermore, the evaluation of H2FLEs for which the weights of criteria are either totally known or unknown is discussed. The BWM approach is used to evaluate the criteria weights for the second circumstance. They are altogether more convincing than the abstract weightdeciding methodology that experts straightforwardly give the weights.

Framework of this study
The overall sequence of this paper is prescribed as follows: some preliminary definitions and existing correlation coefficients are briefly discussed in section 2. The proposed method for H2FLSs with some important properties is discussed in section 3. A brief discussion of BWM is mentioned in section 4. The validity of the work illustrated with the help of example is given in section 5. Comparative analysis with some existing approaches is described in section 6. At long last, the consequence is given at the end.
The continuous LTS S is defined by Xu [49,50], which is an extension of the discrete term set S where S ¼ fu k ju 0 � u k � u g g. If u k 2 S, u k is called the original linguistic term and will be used to assess the alternatives throughout a decision process. If u k = 2 S, then u k is called the virtual linguistic term of S that will appear only in the calculations process. Definition 2 [9] Let x i 2 X, i = 1, 2, � � �, N, be fixed and S = {s t |t = −τ, � � �, −1, 0, 1, � � �, τ} be a LTS. A hesitant fuzzy LTS on X, H S , is in the form of where h S ðx i Þ ¼ fs φ l ðxiÞjs φ l ðxiÞ 2 S; l ¼ 1; 2; :::; Lðx i Þg with L(x i ) being the number of LTS in Herrera and Martínez [12] defined 2-tuple linguistic model to represent the linguistic type of data by 2-tuple (u i , a), where u i 2 S and a 2 [−0.5, 0.5). It means that, a function is transformed into 2-tuple linguistic term and a numerical number.
Definition 3 [12] Let S ¼ fu 0 ; u 1 ; :::; u g g be a LTS and β 2 [0, g] is a numerical value indicating the result of a symbolic aggregation operator. Then, 4 : [0, g] ! S× [-0.5, 0.5) shows a 2-tuple linguistic function that presenting the same information to β is defined as follows: a ¼ b À i; a 2 ½À 0:5; 0:5Þ: ( 4 operator is a one-to-one function and an inverse function is The 2-tuple linguistic model ought to be intelligible and predefined, by which each decision-maker can additionally qualify potential alternatives for specific criteria. However, because of the complication of objects and vagueness of people thinking, the 2-tuple linguistic model can not wholly describe the hesitating circumstances. Therefore, Wei and Liao [23] presented the hesitant 2-tuple fuzzy linguistic data model to deal with the situations in which data portrayed is in linguistic terms, and the decision-maker has some hesitation in choosing its conceivable linguistic interpretations. It can deal with the information viably in decision-making without harming data and more accurate results. Wei and Liao [23] defined the H2FLS as follows.
Definition 4 [23] Let S = {u 0 , u 1 , . . ., u g } be a LTS. A H2FLS on X, H S , is in the form of In our next study, this hesitant 2-tuple fuzzy linguistic data model will be used in our research for the construction of correlation coefficient for H2FLEs.

Mean and variance of H2FLEs
In this section, we established the mean and variance for H2FLEs. The developed mean and variance will further use in the construction of the correlation coefficient for H2FLSs. It can be defined as follows: Definition 5 Let S = {u 0 , u 1 , . . ., u g } be a LTS and h(x) = β i , be a H2FLEs. The mean of h(x) is defined as: where β i = (u i , α i ) and i = 1, 2, . . ., #h. The variance of h(x) is defined as:

Covariance and correlation coefficient of H2FLSs
In the present subsection, covariance and correlation coefficient are discussed. The covariance of two random variables is a numerical measure of the extent to which their values tend to increase or decrease together. As the estimations of covariance depend on the scale of measurement for x and y, and afterwards, it must be standardized before it's utilization like a conventionally practical measure of association, so it is not a reasonable descriptive measure to compute the relationship between two variables, that is the reason, the correlation coefficient is used to measure the relationship between two variables. The measure of correlation is used in many situations for decision-making problems. Therefore, this paper is based on statistics, which measures the correlation coefficient for H2FLSs.
The mean of H2FLS can be defined as: The variance of the H2FLS is defined as: Definition 7 Let S = {u 0 , u 1 , . . ., u g } be a LTS. For two H2FLSs H a S ¼ f< The covariance between H a S and H b S is defined as:

Important properties
ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi VarðH a S Þ p ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi According to Cauchy-Schwartz inequality The most important property of the correlation coefficient lies between -1 and 1, inclusive. We will prove this property for the special case where we have random variables with zero means and unit variances. To prove this property, we standardize the two variables as It is obvious that the expected values or means and variances of the two variables are 0 and 1, respectively. Symbolically, we write them mathematically as: We also know that S ÞÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi

Weighted correlation measures of H2FLSs
It has been observed that arguments about different objects occasionally can be linked with different weights. Hence, while computing the correlation coefficient between H2FLSs, we will consider the influence found from the weights. In the current section, we developed a weighted correlation coefficient in the hesitant 2-tuple fuzzy linguistic environment. Suppose that the weight associated with each objective x i is ω i , where ω i 2 [0, 1] (i = 1, 2, . . ., n) and P n i¼1 o i ¼ 1: 1. The weighted mean of the H2FLS is: The weighted variance of the H2FLS is: 3. The weighted covariance between H a S and H b S is:

Important properties
Now we need to verify that r o ðH a S ; H b S Þ satisfies the following properties.
The most important property of the correlation coefficient lies between -1 and 1, inclusive. We will prove this property for the special case where we have random variables with zero means and unit variances. To prove this property we standardize the two variables as: It is obvious that the expected values or means and variances of the two variables are 0 and 1, respectively. Symbolically, we write We also know that S ÞÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ðco i h bi À cH b S Þðo i h bi À H b S Þ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi P n i¼1 ðco i h bi À cH b S Þ 2 r ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi P n i¼1 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi

Ordered weighted correlation coefficient of H2FLSs
We first rank the alternatives in ordered weighted correlation measures and afterwards weights them regarding their ranking positions. The ordered weighted correlation has been widely examined in [51]. Here, we discuss the ordered weighted correlation measure between H2FLSs. Suppose that the weight associated with each position is γ i , where γ i 2 [0, 1] (i = 1, 2, . . ., n) and P n i¼1 g i ¼ 1: Then, the understudy idea identified with the ordered weighted correlation measure can be presented as: 1. The ordered weighted mean of the H2FLS is: The ordered weighted variance of the H2FLS is: 3. The ordered weighted covariance between H a S and H b S is: The ordered weighted correlation coefficient between H a S and H b S is:

Method to determine the criteria weights
Experts have weighed the measurements in some circumstances, but sometimes it just does not. For the case, we must determine the weight of each criterion where the weight vector is not exact. In a situation where there is no evidence to confirm differences in criteria, we can assume that all criteria are considered equal. Otherwise, we must face a strategy to determine the weight of the criteria. We should handle the strategy to measure the criteria weights in the position where there is no evidence to confirm the difference in patterns. This part applies BWM to solve this issue. Here, we shortly discuss the procedure of BWM that can be utilized to obtain the weights of criteria [52]. Fig 1 shows summary of this procedure.
Step 1. Determine a set of decision criteria: In this progression, the experts recognizes n criteria {C 1 , C 2 , . . ., C n } that are utilized to create a decision.
Step 2. Decide the best (for example most attractive, most significant) and the worst (for example least attractive, least significant) criteria.
Step 3. Decide the preference of the best criteria over all the others, H2FLEs of best criterion to all the other criteria is performed and then transformed to numbers according to (1). The consequent best-to-other vector would be: Q B = (q B1 , q B2 , . . ., q Bn ), where q Bj shows the preference of the best criterion B over criterion j. Obviously q BB = 1 Step 4. Decide the preference of the considerable number of criteria over the worst criteria, utilizing H2FLEs and then transformed to numbers according to (1). The consequent others-to worst vector would be: Q W = (q 1B , q 2B , . . ., q nB ) T , where q jW Shows the preference of the criterion j over the worst criterion W. Obviously q WW = 1 Step 5. Determine the feasible weights ðP � w1 ; P � w2 ; :::; P � wn Þ [52]. To derive the above weights presented in Step 5, we minimize the maximum among the set of {|P wB − q Bj P wj |, |P wj − q jW P wW |}. Hence it can be formulated as follows: min max j fjP wB À q Bj P wj j; jP wj À q jW P wW jg subject to P P wj ¼ 1 ð9Þ P wj ≽ 0; for all j.
Model (9) can be transferred to the following linear programming model: minz L subject to jP wB À q Bj P wj j � z L ; for all j jP wj À q jW P wW j � z L ; for all j P P wj ¼ 1 ð10Þ P wj ≽ 0; for all j By solving problem (10), the optimal weights ðP � w1 ; P � w2 ; :::; P � wn Þ and z L are obtained.

MCGDM problem solving approach based on the proposed correlation coefficient for H2FLSs
In the current subsection, a method of MCGDM for decision-making purposes in the environment of H2FLSs is analyzed. Assume, there is a set of n alternatives Step 1. The configuration of H2FLSs judgment matrices R(t) (t = 1, . . ., K) established in the authoritative perspective on the efficiency of A i (i = 1, . . ., n) on criteria C j (j = 1, . . ., m) using the LTS S. (see [24]).
Step 2. The configuration of H2FLSs as summation findings of matrix R. (see [24]).
Step 3. Check out even if the weight vector of criteria is presented by the decision-makers or not. If present, proceed to Step 5. If not, go to the next step.
Step 4. Set up the weight vector by linear model of BWM [52] (section 4), proceed to the next step.
Step 5. Computation of correlation coefficient between different alternatives with an ideal alternative. An ideal alternative has the maximum membership degree and minimum nonmembership degree.
Step 6. Ranking of the alternatives on the basis of measurement of the correlation coefficient.
The best alternative has the largest correlation coefficient.

Numerical example (with known & unknown criteria weights)
To feature the practicality and viability of the planned decision-making method, two models are utilized from some relevant literature [24,53,54]. Assume an organization needs to pick a provider with the most grounded far quality for improving its administration levels. There involve four suppliers A i = {A 1 , A 2 , A 3 , A 4 }. Three assessment criteria are considered for four goods providers, i.e., quality (C 1 ), supply capacity and flexibility (C 2 ), and price (C 3 ) with ω = (0.504, 0.335, 0.161). Experts give their views based on the LTS, S = {u 0 =nothing, u 1 = very poor, u 2 = poor, u 3 = moderate, u 4 = good, u 5 = very good, u 6 = perfect}.
Step 5. The correlation coefficients of each alternative A i (1 � i � 4) and A � with respect to their given weight are calculated as follows in Table 2.  Step 6. In response to the estimations of correlation coefficients in step 3, the alternatives are ranked. The ranking order is acquired as Consequently, The organization will pick A 4 for their services.
In the given an example, it is seen that the weight vector is presented already by the decision-makers, but sometimes, it is likely that the weight vector is obscure. If so, we can utilize the BWM to determine the weights of the criteria. Here we again go back to step 4 only to show what to do if we lack this knowledge and set up the weight vector.
Step 4. The company recognizes Quality (C 1 ) for best and Price (C 3 ) for the worst criteria. Tables 3 and 4 demonstrate the comparison vectors. Due to the ambiguity and intangibility of decision-making, H2FLSs is used for optimal company selection.
After the simplification of the given data presented in Tables 3 and 4 by using (1), we have the following results in the model (10) for this problem, as follows: min z L subject to Solving this model, we have (stage 5): P � w1 ¼ 0:628; P � w2 ¼ 0:222 and P � w3 ¼ 0:150:  Table 3. Pairwise comparison vector for the best criterion. Criteria Best criterion:

PLOS ONE
Hesitant 2-tuple fuzzy linguistic multi-criteria decision-making method based on correlation measures Step 5. The correlation coefficients of each alternative A i (1 � i � 4) and A � with respect to their weights are calculated as follows: Step 6. As per the estimations of correlation coefficients in Tables 5 & 6, we come to know that the ρ ω (A 4 , A � ) in Table 5 and r g ðA 1 ; A � Þ in Table 6, has the maximum correlation.

Comparative analysis
The area of study being talked about here keep up a comparative study among the existing methods proposed by [24,52,53] with our work. Referring to Table 7 of the above example, the ranking order of the alternatives is as follows: Wang et al. [24] dealt with a similar example utilizing hesitant 2-tuple linguistic Bonferroni mean operator and prioritized weighted hesitant 2-tuple linguistic Bonferroni mean operator, in which the prioritization relationship of specialists is considered. By analyzing the two ranking outcomes in Table 7, the proposed order weighted correlation, and Wang et al. method are similar. However, there exists a difference in positioning consequence of the first two alternatives A 1 and A 4 of the proposed weighted correlation with wang et al. alternatives ranking order. The reason is that the different prioritization relationships of experts considered by Wang et al. can cause changes in overall values. However, we can say that the overall ranking order is stable and adequate. Ge & Wei [53] also coped with a similar example using hesitant 2-tuple weighted averaging operator and hesitant 2-tuple ordered weighted averaging operator, in which the prioritization relationship of specialists is not thought of, and criteria weights are given in advance. Wang et al. [54] likewise proposed a technique to deal with a similar illustration by utilizing hesitant 2-tuple linguistic prioritized weighted averaging operator and hesitant 2-tuple linguistic correlated averaging operator. In [53,54] and our proposed approach, the most feasible option is A 4 , which indicates the validity of the proposed work.
Similarly, it can also be seen that the ranking order obtained from the proposed method using BWM criteria weights is very similar to the ranking order of alternatives obtained in [24,53,54]. From the above comparisons, the obtained ranking orders of the alternatives show the efficiency of the proposed correlation coefficient. The proposed approach can be applied in every field where there is a relationship between two variables. Furthermore, the proposed approach is easy to calculate for H2FLSs compared to other methodologies mentioned in section 1, which are based on different approaches. As a result, we can conclude that our work is strong and capable of solving the MCGDM problem in the context of H2FLSs compared to the existing approaches' computational difficulties.

Conclusion
The correlation measure of H2FLSs investigated in this paper, which lies in the [−1, 1], gives us more information than that the correlation presented by [26]. The presented correlation coefficient in the context of the hesitant 2-tuple fuzzy linguistic environment indicates the degree of relationship between two variables and investigates their linear change in a positive or negative direction. We first explored the ideas of the mean and variance of H2FLEs, and afterwards proposed a new correlation coefficient. Linked to current methodologies, the most important feature of the proposed method identifies the positive and negative associations for H2FLSs. We additionally defined the weighted and ordered weighted correlation coefficients of the H2FLSs. We have proved some desired properties for the covariance and the correlation coefficient. BWM has also been used to determine the weights of criteria and alternatives w.r.t different criteria. The validation of this work is strengthened by the parallel findings of the suggested method with the available ones. The proposed work is appropriate for H2FLSs and easy to calculate, also capable to holds the property of solving MCGDM problems. In the future, we plan to propose the correlation coefficient for intuitionistic H2FLSs and apply them to the MCGDM problems.